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Gamma-ray bursts (GRBs), which have isotropic energy up to 10^^ erg, would be the ideal tool 
to study the properties of early universe: including dark energy, star formation rate, and the metal 
enrichment history of the Universe. We will briefly review the progress on the field of GRB 
cosmology. Meanwhile, X-ray flares, which may have important clues to the central engine, are 
common phenomena in the GRB afterglows. We present statistical results of X-ray flares, i.e., 
energy, duration time and waiting time distributions, and compare the results with solar flares. 
The similarity between the two kinds of flares are found, which may indicates that the physical 
mechanism of GRB X-ray flares is magnetic reconnection. 
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1. Introduction 

The accelerating Universe at the present epoch was discovered from the observation of type 
la supernovae (SNe la) [1,2]. The cosmic accelerated expansion has also been confirmed by other 
observations, such as cosmic microwave background (CMB) [3], and large-scale structure (LSS) 
[4]. But the SNe la can only be detected at low redshifts, i.e., z < 2.0. Gamma-ray bursts (GRBs), 
as the most powerful explosions in the Universe [5, 6], can fill the gap between SNe and CMB. 
GRBs can serves as the complementary tools to measure dark energy and cosmic expansion [7,8, 
9, 10, 11, 12]. Long GRBs can be made “relative standard candles", using luminosity correlations 
that have been found in prompt and afterglow phases [9, 13, 14]. 

X-ray flares wifh shorf rise and decay limes are discovered by Swiff salellife. The isofropic- 
equivalenl energy of X-ray flares is from ~ 10"^^ fo 10^^ erg. The occurrence times of X-ray 
flares range from a few seconds to 10^ seconds after the GRB trigger. Until now, the physical 
origin of X-ray flares has remained mysterious, although some models have been proposed. Here 
we investigate the energy release frequency distribution, duration-time frequency distribution and 
waiting time distribution of GRB X-ray flares [15]. On the other hand, it is well known that solar 
flare is triggered by a magnetic reconnection process. So we will compare these two kind of flares. 

2. Constraints on dark energy 

GRBs have been widely used to constrain dark energy. We use the calibrating method to stan¬ 
dardize 151 long GRBs with ^iso — correlation. Just as using Cepheid variables to standardize 
SNe la, the GRBs can be calibrated with SNe la. We use the latest Union 2.1 data. This method is 
also cosmological model-independent [16]. The full GRB data is separated into two groups. The 
dividing line is the highest redshift in SNe la Union 2.1 data, namely z = 1.414. The low-redshift 
group (z < 1.414) includes 61 GRBs, and the high-redshift group (z > 1.414) contains 90 GRBs. 
The left panel of Figure 1 shows the la and 2a constraints on the plane. Adding the 

GRB data can give much more tighter constraints. 

We also use the nonparametric technique to constrain the uncorrelated equation of state (EOS) 
of dark energy w(z). The EOS is the most important pai'ameters of dark energy. We apply the 
nonparametric method to a joint data set of the latest observations including SNe la, CMB from 
Planck and WMAP polarization, large scale structure, the Hubble parameter measurement and 
GRBs. The right panel of Eigure 1 shows the evolution of uncorrelated EOS of dark energy, which 
is consistent with ACDM at the 2a confidence level [17]. 

3. Star formation rate 

The high-redshift star formation rate (SER) is important in many fields in astrophysics. Direct 
SER measurement is still out of reach by present instruments, particularly at the faint end of the 
galaxy luminosity function. Eong GRBs triggered by the death of massive stars, provide a com¬ 
plementary technique for measuring the SER. In order to test the GRB rate relative to the SER, 
we must choose bursts with high luminosities, because only bright bursts can be seen at low and 
high-redshifts, so we choose the luminosity cut Ljso > 10^^ erg s^* [18] in the redshift bin 0 — 4. 
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Figure 1: Left: la and 2a constraints on and Ha from different data sets. The solid line shows the 
Q.k =0 case. Right: Estimation of the uncorrelated dark energy EOS parameters at different redshift bins 
(wi,W 2 , ...,wio) from SNe Ia+BAO+WMAP9+H(z) data. The open points show the best fit value. The error 
bars are la and 2a confidence levels. The dotted line shows the cosmological constant. Adopted from [17]. 


This removes many low-redshift, low-Liso bursts that could not have been seen at higher redshift. 
The SFR can be estimated as 
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where is the observed GRB counts, and 5 = 0.5 [19]. The derived SFR from GRBs are 

shown as filled circles in the left panel of Figure 2. Error bars correspond to 68% Poisson confi¬ 
dence infervals for fhe binned evenfs. The high-redshiff SFRs obviously decrease wifh increasing 
redshifls, allhough an oscillalion may exisl. We find lhal fhe SFR at z > 4.48 is proportional to 
(1 using minimum method, which is shown as solid line in Figure 2. The right panel 

shows the optical depth Tg due to the scattering between the ionized gas and the CMB photons. The 
WMAP nine-year data gives Tg = 0.089 ±0.014 [20], which is shown as the shaded region. So our 
GRB-inferred SFR can reproduce the CMB optical depth. 


4. Pop III GRBs and metal enrichment 

Absorption lines on the spectra of bright background sources, such as, are main sources of 
information about the chemical properties of high-redshift universe [21]. We simulate the formation 
of first galaxy. When a GRB explodes, its light will go through the metal-polluted region, so 
absorption line will appear in the spectrum. We found that the metals in the first galaxies produced 
by the first supernova explosion are likely to reside in low ionization states (C II, O I, Si II and 
Fe II). The GRB afterglow travel through the polluted environment in the first galaxy. The metals 
will produce the metal absorption lines in the GRB spectrum, as shown in Figure 4. Because the 
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Figure 2: The cosmic star formation history. The dashed line shows their fitting result. The filled circles are 
the SFR derived from GRBs in this work. Adopted from Wang (2014). 
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Figure 3: Possible explosion sites for Pop III bursts. Shown are the hydrogen number density and metallicity 
contours during the assembly of a first galaxy, averaged along the line of sight within the central ~ 100 kpc 
(comoving), at z — 16.4, closer to the virilization of the atomic cooling halo. Now, metals are being re¬ 
assembled into the growing potential well of the first galaxy. The topology of metal enrichment is highly 
inhomogeneous, with pockets of highly enriched material embedded in regions with a largely primordial 
composition. Adopted from [22]. 


topology of metal enrichment could be highly inhomogeneous, so along different lines of sight, 
the metal absorption lines may show a distinct signature. The metal absorption lines, in terms of 
their flux density and corresponding equivalent widths, will be detectable in the near future, with 
the launch of the JWST. From Figure 4, we can distinguish whether the Pop III progenitor dies as 
PISN, or Type II supernova from the observation of Pop III GRBs spectrum [22]. 

5. Self-organized criticality in GRBs 

Thanks to the rapid-response capability and high sensitivity of the Swift satellite, numerous 
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Figure 4: The total GRB spectrum with the metal absorption lines at the reverse shock crossing time (f 0 = 
16.7 X (1 + 16.37)s). Adopted from [22], 




Figure 5: Left: Energy-release frequency distributions, a, The differential energy frequency distribution 
of solar hard X-ray flares, b. The cumulative energy distribution of GRB X-ray flares. The 83 GRB X- 
ray flares are shown as black dots. The black curve gives the cumulative energy distribution N{> E) — 
a + —Fma)?^] with Ue ^ 1.06±0.15 for GRB X-ray flares. Right: Duration-time frequency 

distributions, a. The relation between the occurrence rate and duration time for solar flares, b. The relation 
between the occurrence rate and duration time for GRB X-ray flares shown as green dots. Adopted from 

[15]. 


unforeseen features have been discovered, one of which is that about half of bursts have large, late¬ 
time X-ray flares with short rise and decay times. The unexpected X-ray flares with an isotropic- 
equivalent energy from 10'^^ to 10^^ ergs have been detected for both long and short bursts. We find 
that X-ray flares and solar flares share three statistical properties: power-law frequency distributions 
for energies (left panel of Figure 5), durations (right panel of Figure 5), and waiting times, which 
can be explained by self-organized criticality (SOC) theory [15]. Guidorzi et al. (2015) found that 
the waiting time distributions of gamma-ray pulses and X-ray flares of GRBs also had power-law 
tail extending over four decades, which can be understood by SOC theory [23]. For a SOC system 
[24, 25, 26], owing to some driving force, subsystems will self-organize to a critical state at which 
a small perturbation can trigger an avalanche-like chain reaction of any size within the system. The 
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two types of flares are driven by a magnetic reconnection process. 

6. Conclusions 

GRBs are promising tools for studying the properties of early Universe, including dark energy, 
star formation rate and metal enrichment history. In the future, the French-Chinese satellite Space- 
based multi-band astronomical Variable Objects Monitor (SVOM) and JWST, have been optimized 
to increase the number of GRB and the synergy with the ground-based facilities. This will open a 
new window of the dark ages of the Universe. 
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